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Abstract 

We discuss universality in random matrix theory and in the study of Hamil- 
tonian partial differential equations. We focus on universality of critical behavior 
and we compare results in unitary random matrix ensembles with their counter- 
parts for the Korteweg-de Vries equation, emphasizing the similarities between 
both subjects. 



fH ! 1 Introduction 

> ■ 

It has been observed and conjectured that the critical behavior of solutions to Hamil- 
tonian perturbations of hyperbolic and elliptic systems of partial differential equations 
near points of gradient catastrophe is asymptotically independent of the chosen initial 
data and independent of the chosen equation [221 124j. A classical example of a Hamil- 
^SJ ! tonian perturbation of a hyperbolic equation which exhibits such universal behavior, is 

^ ■ the Korteweg-de Vries (KdV) equation 

ut + Quux + e^Uxxx = 0, e > 0. (1.1) 



If one is interested in the behavior of KdV solutions in the small dispersion limit e — )■ 0, 
it is natural to study first the inviscid Burgers' or Hopf equation ut + 6uux = 0. Given 
smooth initial data u{x, 0) = no(x) decaying at ±oo, the solution of this equation is, for 
t sufficiently small, given by the method of characteristics: we have u{x, t) = u(^(x, t)), 
where ^(x, t) is given as the solution to the equation 

x = 6tuo{0+C (1-2) 

It is easily derived from this implicit form of the solution that the x-derivative of u{x, t) 
blows up at time tc = maxggK(-6M' (g)) ' which is called the time of gradient catastrophe. 
After this time, the Hopf solution u{x, t) ceases to exist in the classical sense. For t 
slightly smaller than the critical time tc the KdV solution starts to oscillate as shown 
numerically in [35]. For t > tc the KdV solution develops a train of rapid oscillations 
of wavelenght of order e. In general, the asymptotics for the KdV solution as e — t- 
can be described in terms of an equilibrium problem, discovered by Lax and Levermore 
[40l [41] . The support of the solution of the equilibrium problem, which depends on x 
and t, consists of a finite or infinite union of intervals |34[ll8j. and the endpoints evolve 
according to the Whitham equations |50ll32j . For t < tc, the support of the equilibrium 
problem consists of one interval and the KdV solution as e — )• is approximated by 
the Hopf solution. For t > tc the support of the equilibrium problem may consists 
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of several intervals and the KdV solution is approximated as e — )• by Riemann 9- 
functions ISH [Ml [20l HI] . 

The (x, t)-plane can thus be divided into different regions labeled by the number of 
intervals in the support of the Lax-Levermore minimization problem. Such regions are 
independent of e and depend only on the initial data. Those regions are separated by 
a collection of breaking curves where the number of intervals in the support changes. 
We will review recently obtained results concerning the asymptotic behavior of KdV 
solutions near curves separating a one-interval region from a two-interval region. The 
two interval region corresponds to the solution of KdV being approximated as e — )• 
by the Jacobi elliptic function, the one interval region corresponds to the solution of 
KdV being approximated by the Hopf solution (jl.2p . 



On the space of n x n Hermitian matrices, one can define unitary invariant proba- 
bility measures of the form 

1 " 

^ exp(-iV Tr V{M))dM, dM = dMu JJ dRe Mijdlm. Mij , (1.3) 

i=i i<j 

where Z„ = Zn{N) is a normalization constant which depends on the integer N and 
V is Su real polynomial of even degree with positive leading coefficient. The eigenvalues 
of random matrices in such a unitary ensemble follow a determinantal point process 
defined by 

(1.4) 

t<] 1=1 

with correlation kernel 

Kn{u,v) = {pniu)pn-l{v) - Pn{v)pn-l{u)) , (1.5) 

U — V Kn 

where pk is the degree k orthonormal polynomial with respect to the weight e~^^ 
defined by 

p,(s)pfc(s)e-^^Wds = 5,k, J, A: e M, 

and Kfc > is the leading coefficient of pk- The average counting measure of the 
eigenvalues has a limit as n = — t- oo. We will denote this limiting mean eigenvalue 
distribution by For a general polynomial external field V of degree 2m, the support 
of consists of a finite union of at most m intervals [16] . If V depends on one or 
more parameters, the measure ^.y will in general also vary with those parameters. 
Critical phenomena occur when the number of intervals in the support of /iy changes. 
A decrease in the number of intervals can be caused essentially by three different events: 

(i) shrinking of an interval, which disappears ultimately, 

(ii) merging of two intervals to a single interval, 

(iii) simultaneous merging of two intervals and shrinking of one of those intervals. 
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Near such transitions, double scaling limits of the correlation kernel are different from 
the usual sine or Airy kernel. At a type (i) transition, the limiting kernel is built out of 
Hermite polynomials |30^ [8l B3 l [T]. at a type (ii) transition the limiting kernel is built 
out of functions related to the Painleve II equation [5l[13], and at a type (iii) transition 
the limiting kernel is related to the Painleve I hierarchy (TIES]. Higher order transitions, 
such as the simultaneous merging and/or shrinking of more than two intervals, can also 
take place but will not be considered here. Rather than on the limiting kernels, we will 
concentrate on the asymptotic behavior of the recurrence coefficients of the orthogonal 
polynomials, defined by the three-term recurrence relation 

SPn{s) = 7„+ip„+i(s) + /3„p„(s) + 7nPn-l(s)- (1-6) 

The recurrence coefficients contain information about the orthogonal polynomials and 
about the partition function Z„ of the determinantal point process (|1.4p ['^ |6l [28] . The 
large n, N asymptotics for the recurrence coefficients show remarkable similarities with 
the asymptotic behavior for KdV solution u{x,t, e) as e — 0. 



2 Phase diagram for the KdV equation 

We assume throughout this section that the (e- independent) initial data uo{x) for the 
KdV equation are real analytic in a neighborhood of the real line, negative, have a 
single local minimum xm for which uo{xm) = —1, and that they decay sufficiently 
rapidly as x — )• oo in a complex neighborhood of the real line. The neighborhood of 
the real line where uq is analytic and where the decay holds should contain a sector 
{| argx| < 5}U{| arg(— x)| < 5}. In addition certain generic conditions have to be valid; 
we refer to [9] for details about those. A simple example of admissible initial data is 
given by uq{x) = — sech^(x). 

2.1 Regular asymptotics for the KdV solution 

Before the time of gradient catastrophe tc = maxggK(-6M' (g)) ' asymptotics for the 
KdV solution u{x, t,e) as e — )• are given by 

u{x,t,e) = u{x,t) + ©(e^), 

where u{x,t) is the solution to the Hopf equation with initial data ^0(2;), i.e. the 
implicit solution uo{^{x,t)) defined by (11.21) . The leading term of the above asymptotic 
expansion was obtained in |40] while the error term was obtained only recently for a 
larger class of equations and initial data in [42] . Such an expansion still holds true after 
the time of gradient catastrophe as long as x is outside the interval where the KdV 
solution develops oscillations. In the oscillatory region, the oscillations for some time 
t > tc can be approximated as e — )• 0, by the elliptic function 

u{x, t, e) = /3i + /32 + /33 + 2a 

+ 2e'^ log ^ (^l^k - 2t(/3i + /32 + ^3) -q];r^+ 0{e). (2.1) 

Here 

+ (...) 

K[s) K[s) I3i - /33 
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where K{s) and E(s) are the complete elhptic integrals of the first and second kind, 
K'{s) = K{VT^), and ^{z; r) is the Jacobi elliptic theta function. In the formula 
(I2.ip the term f3i + f32 + ^3 + 2a is the weak limit of the solution u{x,t,e) of KdV as 
e — 7- and it was derived in the seminal paper |40] . The asymptotic description of the 
oscillations by theta-function was obtained in |49) . A heuristic derivation of formula 
()2.ip without the phase, was first obtained in [37]. The phase q in the argument of 
the Jacobi elliptic theta function (j2.ip was derived in [20]. It depends on /3i,/32,/33 and 
on the initial data and it was observed in |35j that q satisfies a linear over-determined 
system of Euler-Poisson-Darboux type derived in \4:7\ I36|. The negative numbers /3i > 
/32 > /33 depend on x and t and solve the genus one Whitham equations f50]- The 
complete solution of the Whitham equation for the class of initial data considered, was 
derived in ^J. 

At later times, the KdV solution can, depending on the initial data, develop multi- 
phase oscillations which can be described in terms of higher genus Whitham equations 
|32| and in terms of Riemann 6 functions |491 HOl [20] . 

The parameters Pi, (32,^3 can be interpreted in terms of the endpoints of the sup- 
port [0, y/Ps + IJ U [\//52 + 1) + IJ of the minimizer of the Lax-Levermore energy 
functional [iQliniH]. 

A transition from the elliptic asymptotic region to the Hopf region can happen in 
three different ways: 

(i) /?! approaches /32 (shrinking of an interval), 

(ii) (32 approaches /^a (merging of two intervals) , 

(iii) /3i,/?2) and (33 approach each other (simultaneous shrinking and merging of inter- 
vals). 

The transitions (i), (ii) and (iii) will lead to an asymptotic description of the KdV 
solution which is similar to the asymptotic description for the recurrence coefficients 
of orthogonal polynomials when the number of intervals in the support of the limiting 
mean eigenvalue density of random matrix ensembles changes. A transition of type (iii) 
takes place at the point of gradient catastrophe. In the {x, t) plane after the time of 
gradient catastrophe, the oscillations asympotitcally develop in a V-shape region that 
does not depend on e see Figured) At the left boundary (the leading edge), a transition 
of type (ii) takes place, and at the right boundary (the trailing edge) we have a type 
(i) transition. Given t sufficiently short after the time of gradient catastrophe tc, the 
leading edge x~{t) is characterized by the system of equations 

x-{t) = 6tu{t) + fL{uit)), (2.3) 
6t + e{v{t);uit)) =0, (2.4) 
dj{v{t);u{t)) =0, (2.5) 

where u{t) > v{t), fhiu) is the inverse of the decreasing part of no(x), and 9 is given 
by 

This corresponds to the confluent case where the elliptic solution (12. ip degenerates 
formally to linear oscillations, namely (32 = (33 = v and (3i = u. The trailing edge on 
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Figure 1: Sketch of the phase diagram for the equilibrium problem associated to the 
KdV equation. Outside the cusp-shaped region, the support of the Lax-Levermore 
minimizer consists of one interval; inside of two intervals. At the cusp point, we have a 
type (iii) transition, at the left breaking curve one of type (ii), and at the right breaking 
curve, one of type (i). 



the other hand is characterized by 



x+{t) = 6tu{t) + fL{u{t)), 
6t + e{v{t);u{t)) = 0, 

v{t) 



u{t) 



{6t + 9{X; u{t)))y/X-u{t)d\ = 0, 



(2.7) 
(2.8) 

(2.9) 



with u{t) < v{t), and 9{\;u) defined in (|2.6p . In this case we have Pi = (32 = v and 
/33 = u. In this case the solution (12. ip degenerates formally to a soliton. 



2.2 Critical asymptotics for the KdV solution 
2.2.1 Point of gradient catastrophe 

Near the first break-up time, the KdV solution starts developing oscillations for small 
e. These oscillations are modeled by a Painleve transcendent U{X,T), defined as the 
unique real smooth solution to the fourth order ODE 



X = TU 



1 



1 



6 24 
with asymptotic behavior given by 

C/(X,T) = T(6|X|)i/3 



7:U' + -rm + 2UUxx) + WTT^UxXXX 
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^-Q^I^T\X\-'^'^ + 0{\X\-^), 
3 



(2.10) 



asX^iboo, (2.11) 



for each fixed T G M. The existence of a pole free solution of (j2.10p with asymptotic 
conditions (j2.1ip was conjectured in p2]) and proved in [14j. Let us denote tc for 
the time of gradient catastrophe, Xc for the point where the x-derivative of the Hopf 
solution blows up, and Uc = u{xc,tc)- We take a double scaling limit where we let 
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e — 7- and at the same time we let x — )• Xc and t — )• in such a way that, for fixed 
X,T eR, 

X — Xc — GuJt — tr) , 6(t — tc) ^ , , 

m ^ ^ = X, m , \ ' = T, (2.12) 

where 

In this double scaling limit the solution u{x, t, e) of the KdV equation (jl.ip has the 
following expansion, 

nix, t,e)=n.+ ( ^) U ( ^ " " " '-^] + O (e^^A . (2.13) 

The idea that the solution of KdV near the point of gradient catastrophe can be ap- 
proximated by the solution of (|2.10p appeared first in |441 H5] and in a more general 
setting in [22], and was confirmed rigorously in [9j. In [12] the correction term of order 
e'^/'^ was determined. 

2.2.2 Leading edge 

Near the leading edge, the onset of the oscillations is described by the Hastings-McLeod 
solution to the Painleve II equation 

q"{s) = sq+2q^{s). (2.14) 

The Hastings-McLeod solution is characterized by the asymptotics 

q{s) = y^_s/2(l-Fo(l)), ass^-oo, (2.15) 

q{s) = Ai{s){l + o{l)), ass^+oo, (2.16) 

where Ai(s) is the Airy function. The leading edge x~{t) is, for t sufficiently short after 
tc, determined by the system of equations (|2.3p - (|2.5p . Let us consider a double scaling 
limit where we let e — t- and at the same time we let x — t- x~{t) in such a way that 

lim^^|3^=X GR, (2.17) 

for t > tc fixed. In this double scaling limit, the solution u{x, t, e) of the KdV equation 
with initial data uq has the following asymptotic expansion, 

u{x, t,e) = u- ^q [six, t, e)] cos (^^^ + 0(e§ ), (2.18) 



where 



e{x,t) = 2V^r^{x-x-) + 2 f ifLiO + ^t)VC-vdC, (2.19) 



and 



92 



c = -Vu - v-^9{v; u) > 0, s{x,t,e) = — , — ^ ry^r, (2.20) 
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with 9 defined by (j2.6|) . and q is the Hastings-McLeod solution to the Painleve II 
equation. Here x~ and v < u (each of them depending on t) solve the system ()2.3p - 
()2.5p . The above result was proved in [10], confirming numerical results in [35]. In [lOj . 

2 

an explicit formula for the correction term of order e s was obtained as well. We remark 
that a connection between leading edge asymptotics and the Painleve II equation also 
appeared in [38] . 



2.2.3 Trailing edge 

The trailing edge x~^{t) of the oscillatory interval (i.e. the right edge of the cusp-shaped 
region in Figure 2) is determined by the equations ()2.7p - (|2.9p . As e — t- 0, we have, for 
fixed y and t, 

x+ + — ^^^y, t, e ) =n + 2(v-n)ysech2(Xfc) + 0(eln2e), (2.21) 

where 

Xk = ^{^-y + k)liie- ln{V2^hk) " + ^) l^T, 

2f , (2-22) 

hk= 1 7 = 4(7; - n) 4 y/-d^0{v; u), 

vr* V A:! 

and ^ is given by ()2.6p |llj . It should be noted in this perspective that the KdV 
equation admits soliton solutions of the form a sech^(6x — ct). This means that the last 
oscillations of the KdV solution resemble, at the local scale, solitons. 



3 Phase diagram for unitary random matrix ensembles 
3.1 Equilibrium problem 

In unitary random matrix ensembles of the form (jl.3p . the limiting mean eigenvalue 
density is characterized as the equilibrium measure minimizing the logarithmic energy 

Ivif^) = log j-^^dfi{s)dfi{y) + I V{s)df,{s), (3.1) 

among all probability measures on M. For a polynomial external field of degree 2m, the 
equilibrium measure is supported on a union Sy of at most m disjoint intervals. Its 
density can be written in the form [IB] 



i'vis) = n V (^i - ^)(« - «i) ^(^)' s ^ ^j=i[aj:bj], k < m, (3.2) 

where h is a polynomial of degree at most 2(m — k). The equilibrium measure is 
characterized by the variational conditions 

2 j log|s - y\dfi{y) - Vis) = iv, s€ U^j^^[aj,bj], (3.3) 

2 log\s-y\dn{y)-V{s)<iv, s G M. (3.4) 



The external field V is called /c-cut regular if h{s) in (13. 2p is strictly positive on 
UJ^-^^[aj, bj] and if ()3.4p is strict for s G M \ Uj^^laj, bj]. In other words, it is singular if 
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(i) equality in (|3.4|) holds at a point s* € M \ L)j^-^[aj,bj], 

(ii) h{s*) = with s* G u'^^-^{aj,bj). 
(in) h{s*) = with s* = aj or s* = bj. 

3.2 Example: quartic external field 

Let us now study a two-parameter family of quartic external fields 

T4j(s) = e^ (1-t) ht ^ \--s 

^ ^ ^21 20 15 5 5 



(3.5) 



For t = 0, we have V^.,o('5) = e^^, which means that the random matrix ensemble is a 
rescaled Gaussian Unitary Ensemble. The equilibrium measure fix.o is then given by 



d/ix o{s) = — \/4e-^-s2ds, s e [-2e-^/2^ 2e-^/2]_ (3 g) 

' zvr 

It can indeed be verified directly that this measure satisfies the variational conditions 
(!3.3p - (|3.4p . For x = and < t < 1, one can verify that 



"'/"o,t(g) = ^^^^^^^2) Vs^-ms-2f + -f^)ds, sG[-2,2], 7 = ^^-5.(3.7) 

This shows that Vo,i has a singular point of type (iii) at s = 2. On the line t = 9, 
Vx^g{s) is symmetric around s* = |. The external field is one-cut regular for x < x* =: 
— log and presumably two-cut for x > x*. For x < x*, the equilibrium measure is 
given by 



"^^"'^^""^ " 7r62(62^+4g) V^^ " ^ + + ^ " s)((s-^)2+C)ds, s G 

(3.8) 

where 



b=\l^ + ^V5e-V27e- + 245e2- C = ^(80-96V) (3.9) 
At x = X*, the equilibrium measure is given by 



dl^x',9{s) = :^^J{s-^ + b)i^ + b-s)is-^fds, ^^[^-^4^^^' ^=^^' 

which means that there is a type (ii) singular point at s* = 4/3. 

For t fixed and x sufficiently large and positive, it follows from results in |39] that 
the number of intervals is equal to the number of global minima of Vx^t, which is one 
for i < 9 and two for t = 9. For t fixed and x sufficiently large negative, one can show 
that the equilibrium measure is supported on a single interval. Also, for any t, when x 
decreases, the support of the equilibrium measure increases. This suggests that there 
are, as shown in Figure [21 two curves in the (x, t)-plane where Vx^t is singular: one 
connecting (0,1) with (2;*, 9) where a singular point of type (ii) is present, and one 
connecting (0, 1) with (+00, 9) where a singular point of type (i) occurs. 

Remark 3.1 In [2j, orthogonal polynomials with respect to complex weights of the 
form e-"-^(^) were considered, with V quartic symmetric with complex-valued leading 
coefficient. This lead to a phase diagram which shows certain similarities with ours, 
but also with breaking curves of a different nature. 
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Figure 2: Sketch of the phase diagram for the equihbrium measure in external field 
Vx^t- The one-cut region and the two-cut region are separated by two curves, at the left 
curve a type (ii) singular point is present, at the right curve a type (i) singular point, 
and at the intersection point (0, 1) there is a type (iii) singular point. 



3.3 Regular asymptotics 

If y is a one-cut regular external field, the leading term of the asymptotics for the 
recurrence coefficients depends in a very simple way on the endpoints a and h: we have 

m 

7n = ^-j^ + ©(n"^), asn-;>oo, (3.10) 

/3„ = + C'(n"2), asn^oo. (3.11) 

If y is a two-cut regular external field, the leading order term in the asymptotic ex- 
pansion for the recurrence coefficients is still determined by the endpoints ai, 6i, a2, 62; 
but the dependence is somewhat more complicated, and the leading term is oscillating 
with n. An explicit formula for the leading order asymptotics was given and proved in 
[17J for A;-cut regular external fields V ^ with k arbitrary. We will not give details about 
those asymptotics, but we note that the expansion is of a similar nature as (12. ip in the 
two-cut case. 

3.4 Critical asymptotics 

We will now describe the critical asymptotics for the recurrence coefficients 7„(x,t) and 
/3n{x,t) of the orthogonal polynomials with respect to the weig ht e-"^->*. It should be 
noted that critical asymptotics near type (ii) and type (iii) singular points are known 
for more general deformations of external fields Vx^t than only the one defined by (j3.5p . 

3.4.1 Singular interior points 

Assume that T4*,t*(s) is a singular external field with a singular point s* of type (ii) (a 
singular interior point), and with support [a, b] of the equilibrium measure. Asymptotics 
for the recurrence coefficients were obtained in [5j for quartic symmetric V and in |13] 
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for real analytic V. Let us specialize the results to our example where ^ is given by 
()3.5p . Since Vx*^t* is quartic, this implies that ipx*,t* has the form 

^x',t* {s) = C^/{s-a){h-s){s - s*f. (3.12) 

Then as n — t- oo simultaneously with x — ?■ x* such that x — x* = 0{n~'^^^), we have the 
asymptotic expansions [13] 



7„(x, t*) = ^—^ - {sx,n) cos(2™(x))n-i/3 + 0{n-^/^), (3.13) 

/5n(x, t*) = ^ + -q {sx,n) sin(27rna;(x) + 0)n-V3 + 0{n'^'''), (3.14) 

where 

1 



n 



cy^ (s* — a){b — s*) ' 
and where c, 0, and are given by 

1/3 



T,C^{s*-a){h-s*) \ . b + a 

c = , = arcsm , 

\ 4 / b — a 

io{x) = / ijjx*,t*{s)ds + 0{n^'^^^), as n ^ cxD. 



An exact formula for to can be given in terms of a modified equilibrium problem. When 
x approaches x* , we observe that the recurrence coefficients develop oscillations. The 
envelope of the oscillations is described by the Hastings-McLeod solution q. One should 
compare formulas ()3.13p - ()3.14p with ()2.13p and note that the scalings correspond after 
identifying e with 1/n. 



3.4.2 Singular edge points 

Asymptotics for the recurrence coefficients for general one-cut external fields V with 
a singular endpoint were obtained in [15j. Let Vq be an external field such that the 
equilibrium measure is supported on [a, b] and such that the density "00 behaves like 
i^o{s) ~ c{b — s)^/^ as s — )■ 6, c / 0. Double scaling asymptotics were obtained for 
external fields of the form Vq + SVi + TV2 with real S,T — )• 0, where Vi is arbitrary 
and V2 satisfies the condition 

We can write Vx^t in the form 

Vx,t{s) = Vo,i{s) + (e" - l)Fo,i(s) + e"(t - l)(^o,i(^) " ^o,o(s)). (3.15) 

Since 
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we can apply the results of [15]. In the double scaling limit where n — )• oo and simul- 
taneously j; — 7- 0, t — >• 1 in such a way that limn^/'^(e^ — 1) and lim n^/'^'e^ (i — 1) exist, 
we have 

7„(x, t) = l + ^ C/(cin6/7(e- - 1), csn^/'e^Ct - l))n-2/7 + ©(n'^A), (3.I6) 
2c 

t) = - [/(cin6/7(e^ - 1), C2n^/^e^(t - l))n~^/'^ + ©(n"^/^). (3.17) 
c 

The constants c, ci, C2 are given by 
c = 6^/^ > 0, 



C2 



1 f 2 + u 



47ric3/2 7^ Y (2 - 7X)3 
where 7 is a counterclockwise oriented contour encircling [—2,2]. 

Remark 3.2 Applying the results from [15| directly, one has an error term 0{n~^^'^) 
in (j3.16p and (|3.17p . but going through the calculations, it can be verified that the error 
term is actually 0(n~^/''). The analogy between (|3.16p - (|3.17p and (|2.13p is obvious. 

3.4.3 Singular exterior points 

Asymptotics for the recurrence coefficients in the vicinity of a singular exterior point 
have not appeared in the literature to the best of our knowledge. Asymptotics for 
orthogonal polynomials associated to an external field V with a singular exterior point 
and for the correlation kernel (jl.Sp have been studied in [HI [H |33] using the Riemann- 
Hilbert approach. We are convinced that the same analysis can be used, with some 
additional effort, to compute asymptotics for the recurrence coefficients. If is an 

external field with a singular exterior point, the analogy with the KdV asymptotics 
suggests asymptotic expansions of the form 



7n{x*-y—,t) = ^i^!l^L^(^^ + C1 Vsech2(Xfc) + 0(n-iln2n) (3.18) 
con 4 

Mx* - y—,t) = + + g sech^X,) + 0{n-' W n), (3.19) 



Cnn 

^ A;=0 



as n — )■ 00, where 

Xk = -C2{y,k)lnn + c^{k). 



4 The problem of matching 

Asymptotic expansions for KdV solutions are known in the regular regions and in criti- 
cal regions, but we do not have uniform asymptotics for u{x, t, e) in x and t. Indeed, the 
critical asymptotics are only valid in shrinking neighborhoods of the breaking curves: 
a neighborhood of size near the leading edge, a neighborhood of size 0{elne) 

near the trailing edge, and a neighborhood of size 0{e^^'^) at the point of gradient 
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catastrophe. On the other hand, the regular asymptotics are only proved to hold uni- 
formly for X and t at a fixed distance away from the breaking curves. However one 
can see easily that (j2.13p and (|2.2ip match formally with the regular asymptotics for 
X close to the breaking curves but outside the cusp-shaped region. Indeed for (j2.13p 
this follows from the decay of the Hastings-McLeod solution q at+oo. When x is close 
to the boundary but inside the cusp-shaped region, the situation is more complicated. 
One can hope that the regular asymptotics can be improved in such a way that they 
hold also when x, t approach a breaking curve sufficiently slowly when e tends to 0, and 
that the critical asymptotics can be improved to hold in a slightly bigger neighborhood 
of the breaking curves. It would be of interest to see if such an approach could provide 
uniform asymptotics for the KdV solution as e — )• 0. 

The problem of obtaining uniform asymptotics in x and t for the recurrence coeffi- 
cients 7„(x,t) and (3n{x,t) may seem an artificial one at first sight, since one is often 
interested in a random matrix with a fixed external field V instead of letting V vary. 
However, it becomes more relevant when studying the partition function 



where kj is the leading coefficient of the normalized orthogonal polynomial pj with 
respect to the weight e~"^. A consequence of this formula is that, if one lets V vary 
with a parameter r in a convenient way, it is possible to derive various identities for 
T-derivatives of In Zn in terms of the recurrence coefficients 7^ (r) and /3fc (r) for k large 
[6l [28]. A possible strategy to obtain asymptotics for the partition function, is to let 

2 

the r-dependence be such that V interpolates between the Gaussian V{z;to) = ^ 
and V{z;ti) = V{z). Integrating the differential identity then requires asymptotics 
for the Gaussian partition function (which are known) and uniform asymptotics for 
the recurrence coefficients jn{T) and l3n{T) over the whole range [ro,Ti]. Depending 
on the chosen deformation, this could require uniform asymptotics for the recurrence 
coefficients near a singular point of type (i), (ii), or (iii). The results presented in 
the previous section do not provide sufficiently detailed asymptotics for the recurrence 
coefficients: they are not uniform near the breaking curves. For example near a critical 
point of type (ii), formulas ()3.13p - (|3.14p are only valid for x — x* = 0{n~'^^^) asn —?■ 00, 
whereas the asymptotic formula in the two-cut region is valid only at a fixed distance 
away from a critical point. 



It is well-known that recurrence coefficients for orthogonal polynomials follow the time 
flows of the Toda hierarchy. In this section, following ^23, i25j we will formally derive 
the KdV equation as a scaling limit of the continuum limit of the Toda lattice. This 
gives a heuristic argument why asymptotics for KdV and the recurrence coefficients 
show similarities. 




n 



It is well-known that 



5 The Toda lattice and KdV 
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n G 



(5.1) 



The Toda lattice is a Hamiltonian system described by the equations 

dun dvn 

The Toda lattice is a prototypical example of a completely integrable system |31j . Let 

2d 

y{i) = Vo{i) + Y.tji\ t2d>0, (5.2) 

J=l 

where Vo(.^) is a fixed polynomial of even degree with positive leading coefficient, and 
let pj be the orthogonal polynomials defined by 



1 



(5.3) 



where e = — is a small positive parameter. As mentioned before, the polynomials 

Pn{C) satisfy a three term recurrence relation of the form (|1.6p . 

The recurrence coefficients 7^ and in ()1.6p evolve with respect to the times tj. 
defined in (|5.2|) according to the equations [221 ED 133 E] 



- Y [[Q ]n-l,n-l 



dtk 



ln[Q^]n,n-l 



n+l,n5 



(5.4) 
(5.5) 



where [Q ]n,m denotes the n, m-th element of the matrix Q and Q is the tridiagonal 
matrix 



Q 



//3o 


71 











• A 


71 


/5i 


72 













72 


/32 


73 













73 


/33 


74 













74 
















■ •/ 



(5.6) 



The equations (|5.4p - ()5.5p are the Toda lattice hierarchy in the Flaschka variables |31j . 
In particular the first fiow of the hierarchy takes the form 

- KPn-l - Pn) , 

(5.7) 



dPn 2 2 

7n - 7n+l 



These equations correspond to the Toda lattice ()5.ip by identifying ti = t , f3n = —Vn 
and 



Un = log7„. 



(5.^ 



In addition to the Toda equations, the recurrence coefficients for the orthogonal poly- 
nomials satisfy a constraint that is given by the discrete string equation which takes 
the form 1331 



ln[y' {Q)]n,n-l = ne, 
[V'{Q)]n,n = 0. 



(5.9) 
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For example, choosing Vo(C) = ^ one obtains 



= 0) = 0, 7n(t = 0) 



ne, 



(ti, t2, . . . , t2d)- 



(5.10) 



To obtain the continuum hmit of the Toda lattice, let us assume that u{x) and 
v{x) are smooth functions that interpolate the sequences Un,Vn in the following way: 
u{en) = Ufi and v{€n) = f„ for some small e > 0, n > 0, x = en. Then the Toda lattice 
()5.ip reduces to an evolutionary PDE of the form \27\ [T9] 



1 1 

= - [v{x) - v{x - e)] = Vx- -evxx + O(e^) 



1 

vt = - 
e 



^u{x+e) _ u{x) 



1 « 



+ 0(6^). 



(5.11) 



In order to write the continuum limit of the Toda lattice in a canonical Hamiltonian 
form, following Dubrovin-Zhang [26j, we introduce w{x) by 



w(x) = edx[l - e'^"'-']' u{x) = w + -w^ + —w^x + 



(5.12) 



In the coordinates v, w the continuum limit of the Toda lattice equations takes the 
form 



Wt = Vx 



Vt 



Wx + —{2Wxxx + ^WxWxx + wl) 



+ 0(e^ 



with the corresponding Hamiltonian given hy H = J 



— + e^wl + 

2 24 



(5.13) 



dx and 



Poisson bracket {v{x),w{y)} = 5'{x — y) where 5{x) is the Dirac 5 function. We remark 
that in these coordinates the continuum limit of the Toda equation contains only even 
terms in e. For e = 0, (|5.13p reduces to 



wt = Vx, vt = e Wx 



(5.14) 



The solution of equations (j5.14p can be obtained by the method of characteristics. The 
initial data relevant to us should satisfy the continuum limit of the string equation (15. 9p 
for t = 0. The Riemann invariants of ()5.14p are 

r± = V ziz 2e'2 

so that (j5.14p takes the form 



da 

— r± + A±— r± = 0, X± 



=Fe: 



=F- 



r+ 



The generic solution of (j5.14p can be written in the form [48\ (50] 

X = A±t + /±(r+,r__), 
where f± (r+ , r_ ) are two functions that satisfy the equations 

/+ - /- d 



dr 



d _dx+u-f. 



5r_ A+ — A- 



2(r+ — r_) dr- 



-f- 



(5.15) 



(5.16) 
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From the above relations one can conclude that there exists a function / = f{r^,r_ 
so that 



fi 



df_ 

dr+ 



The explicit dependence of / for a certain class of initial data can be found in [19] . To 
obtain / in the random matrix case we impose that the equations (jS.lSp are consistent 
with the continuum limit of the discrete string equation (I5.9P for ti = t > and tj = 



for j > 1- At the leading order in e the string equation (j5.9p in the Riemann invariants 
r± = — /3ib27 gives after straightforward but long calculations, the following expression 
for the function f{r+,r-): 



(5.17) 



Remark 5.1 The equations (j5.15p with / given in (j5.17p . coincide with the equations 
that define the support of the equilibrium measure for the variational problem 



inf 



JrJr \C-m X Jr 



in the case where the equilibrium measure is supported on one interval. The Riemann 
invariants r_|_ and r_ can thus be interpreted as the end-points of the support of the 
equilibrium measure. 

In what follows, we are going to show that the solution of the equation (15.13P in the 
vicinity of a singular point of type (iii) reduces to the KdV equation, in agreement with 
|23j . First we consider the solution of the hodograph equation (]5.15p near a singular 
point of type (iii); namely, let {xc, tc) be a point of gradient catastrophe for the Riemann 
invariant r+, which means that dxV^ goes to infinity at the critical point {xc,tc)- We 
define r±{xc,tc) = r'^- Such a critical point is characterized by the conditions 

^+,+tc + /+,+ = 0, X^^^^tc + /+,++ = 0, 
and the critical point is generic if 

\l,+++t, + / 0, Ai,_t, + /f / 0, 

d 

where in the above formulas we used the notation Al „ = — — A_(r i = r'j_,r_ = r'l) 

and consistently for the other terms. 

Expanding in power series (I5.15P near (xc, tc) and using (I5.16P after the rescalings 



x^ 
f- 



k ^/^(x - Xc - Ai(t - tc)), x+ = k ^(x - Xc - A+(t - tc)) 



(5.18) 



one obtains, letting /c — 0, 
x_ = Clf_ 

X+ = C2X_f+ + csfi. 



(5.19) 



15 



where 



ci = (/V + Ai,_tc), C2 = ^^^, C3 = ^(A;,+++tc + /;,+++). (5.20) 

Wc observe that ()5.19p describes a Withney singularity in the neighbourhhood of (0, 0) 
|23j. Performing the same rescahngs (|5.18p to the equations (j5.14p and letting A; — )• 
one obtains 

9f_ (9f_|_ _ 5f_|_ 

a — = 0' ^ ^'^2r+^ — = 0, 

ax_|_ OX- ax+ 

with C2 as in (j5.20p . Clearly the equations (|5.19p are a solution of the above equations 
with singularity in (x+ = 0, x_ = 0) and at f-|- = 0. The next step is to perform the 

7 

rescaling (|5.18p to the equation (|5.13p and letting e — >• kee. One obtains in the limit 

dxC^' ^ 192(A-_ - A^) = 96 

, (5.21) 
3r+ _ 9r+ 2 ^ - 

^ ^C2?'+T^ hC4e ^-^r+ = 0. 

OX- ox^ Ox^ 

The first of the above equations has been obtained after integration with respect to x^ 
using (|5.19p . The second one is the KdV equation for r+ with time variable x_ and 
space variable x+. Such derivation has been obtained in a more general setting in [25] . 
On the formal level, the above calculations explain why the asymptotic behavior of 
the solution of the continuum limit of Toda lattice and in particular of the recurrence 
coefficients of orthogonal polynomials near the point of gradient catastrophe is of a 
similar nature as the KdV case. However, a rigorous proof of the generic behavior of 
the solution of the continuum limit of Toda lattice near the point of gradient catastrophe 
cannot be derived from the KdV case but a separate proof is needed. 
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